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ABSTRACT 

Because of the possibility of adverse interaction between the control 
system and the structural dynamics of large, flexible spacecraft, great 
care must be exercised to ensure stability and system performance. 
Because of the high cost of insertion of mass into low earth orbit, it is 
prudent to optimize the roles of structure and control systems 
simultaneously. Because of the difficulty and the computational burden in 
modeling and analyzing the control/structure system dynamics, the total 
problem is often split and treated iteratively. The awkwardness and 
inaccuracy of this approach can lead to poor designs. 

It would aid design if the control/structure system dynamics could be 
represented in a single system of equations. Heretofore this has not been 
possible. With the availability of the software PDEMOD it is now possible 
to optimize structure and control systems simultaneously. The 
distributed parameter modeling approach enables embedding the control 
system dynamics into the same equations for the structural dynamics 
model. By doing this the current difficulties involved in model order 
reduction are avoided. 

The NASA Mini-MAST truss is used as an example for studying integrated 
control/structure design. Comparisons are made for (1) structure 
without control, (2) controls only, and for optimal combinations of 
structure and control. Both proof-mass and torque-wheel actuators are 
considered. The results give insight with regard to the essential factors 
in trading structure and control for space applications. 


INTRODUCTION 

Certain future space missions will be performed by large, flexible 
spacecraft. Because of the high cost of insertion of mass into low earth 
orbit, it is prudent to optimize the design of both the structure and 
control systems. The current practice is to create a finite element model 
of the structure and to use a reduced order modal model for control 
syn thesis. Unfortunately, this disjoint procedure is inaccurate and is an 
impediment to integrated design. There is a need for a dynamics model 
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which includes both structural and control dynamics which is 
parameterized in terms of the design variables. 


Distributed parameter modeling has been shown to be quite accurate for 
modelling (reference 1.) the dynamics of the first six modes of portal 
frames. The root-mean- square error in frequency was about 7/10 or one 
percent for the first six modes of two experimental configurations. The 
difficulty in modeling complex structures has been an obstacle to the use 
of distributed parameter modeling. The development of finite element 
software has resulted in the wide-spread practice of modeling dexib e 
structures with finite element models. The availability of modeling 
software such as DISTEL (reference 2.) and BUNVIS-RG (reference 3 J 
offers the alternative of modeling complex configurations with distributed 
parameter models. 


Another advantage of distributed parameter modeling is that control and 
sensor dynamics can be incorporated into the equations of motion ot the 
structural dynamics. Again, the burden of assembling the necessary 
equations for complex configurations has become routine because of the 
capabilities of the software PDEMOD (reference 4). It is now possible to 
optimize structure and control systems simultaneously (references 5. and 
6.) for complex spacecraft because of the incorporation of their dynamics 
into a single system of equations. The distributed parameter modeling 
embeds the control system dynamics into the same equations for the 
structural dynamics model so that model order reduction is not 
necessary. The resulting structures/controls model is particularly well 
suited for integrated design. 


In this paper only preliminary comparisons are made for the Mini-MAST 
truss (references 7. and 8.) for (1) structure without control, (2) controls 
only, and for optimal combinations of structure and control. Both proot- 
mass and torque-wheel actuators are considered. The results give insight 
into the essential factors in trading structure and control for space 
applications. Subsequent, more in-depth study will consider the 
transient dynamics aspects using PDEMOD and will use the optimization 
techniques of Fogel and Holland (references 9. and 10.). 


DISCUSSION 

The software PDEMOD enables the generation of models of complex 
structural configuration which include the dynamics of the control 
system as well as the structural dynamics. This is done using partial 
differential equations to describe the dynamics of flexible beam elements 
which together with rigid body elements form a connected network oi 
components of the structure. The coefficients of the sinusoidal and 
hyperbolic functions for each flexible element give the mode shapes. The 
sensed motion and control forces and moment are expressed in terms oi 


the same parameters and the influence of control on the configuration is 
represented by terms added to the equations for the structural dynamics. 


Structural Dynamics 

First, the structural configuration is viewed as an assembly of flexible and 
rigid elements. It is then necessary to indicate the connectivity of the 
elements. This is done by giving the identification of the rigid bodies at 
enC * eac ^ flexible beam element and the related points of 
attachment. The alignments of the flexible beam axes must also be given 
The axis for each rigid body is that of a particular, attached beam. 

The next input needed is the stiffness (EI X , EI y , EA Z , EI y ) and mass (m/L, 
I y /L] characteristics and the length of each flexible beam element. The 
mass and inertia of each rigid body is needed to complete the 
information required for the structural dynamics model. 


Control Dynamics 

The addition of feedback control does not increase the order of the 
system matrix unless the applied force or moment is applied to the 
interior of a flexible beam element. In such a case it is necessary to add a 
node or rigid body with negligible mass at that point. In other cases it is 
only necessary to add terms to the matrix elements which already exist. 
The dynamics of the sensors and actuators are inserted as transfer 
functions multiplying the additional terms. The additional terms are 
located in the rows which correspond to the body to which the control 
iorce and/or moment is applied, and in the columns which correspond to 
the beam elements which contain the location of the sensed motion. 


Optimization 

In order to perform optimization for parameter estimation or for criteria 
involving structural dynamics, sensitivity functions are usually required. 
T^ie sensitivity functions relate the change in the criterion to changes in 
the parameters involved in the optimization. Although it is possible to 
express in closed form these derivatives it is most unwieldy to evaluate 
the expressions. This is because derivatives must be taken of the 
determinant of the system matrix which can be quite large. It is more 
practical to approximate the derivative numerically by the ratio of the 
change in the criterion to the change in the parameter. This was the 
approach used in reference 8. for a parameter estimation application 
Changes in the modal frequencies for each of the model parameters was 
generated numerically. 
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Although parameter estimation is an example of optimization, it is the 
optimal, integrated control/structural design that we now turn our 
attention. The selection of the design criterion and the corresponding 
conditions or constraints are of paramount importance. An ill-posed 
problem can easily result in nonsensical results which bear no relation to 
the actual design problem. Perhaps the most suitable design criterion for 
the integrated control/structure problem would be the life-cycle-cost of 
the entire system, subject to a set of system performance and structural 
specifications. In many cases it suffices to consider only the total 
structure and control system mass. 

There are alternatives to optimization schemes which require the 
sensitivity functions mentioned earlier. Genetic algorithms (GAs), as 
introduced by Holland (reference 10), are one form of directed random 
search. The form of direction is based on Darwin's theory of the "survival 
of the fittest". In GA's a finite number of candidate solutions or designs 
are randomly (or heuristically) generated to create an initial population. 
This initial population is then allowed to evolve over generations to 
produce new, and hopefully better designs. The basic conjecture behind 
GA’s is that evolution is the best compromise between determinism and 
pure chance. GA’s have the capability to solve continuous, discrete, and 
mixed optimization problems. 

There are four main operations in a basic GA: evaluation, selection, 

crossover, and mutation. Evaluation is the process of assigning a fitness 
measure to each member of the current population. The fitness measure 
is typically chosen to be related to the objective function which is to be 
maximized. No gradient or auxiliary information is used. Therefore, GA's 
are more likely to converge to a global maximum than a hill climbing 
algorithm, although no algorithm can guarantee convergence to the global 
maximum. Selection is the operation of choosing members of the 
current population to be parents for producing the next generation. 
Selection is weighted towards the more fit members of the population. 
Therefore, designs which are better as viewed from the fitness function, 
and therefore the objective function, are more likely to be chosen as 
parents. Crossover is the process in which design information is 
transferred to the prodigy from the parents. Mutation is a low probability 
random operation which slightly perturbs the design represented by the 
prodigy. The mutation operation is used to retain design information 
over the entire domain of the design space during the evolutionary 
process. 


INTEGRATED DESIGN PROBLEM 

The integrated control/structural design problem to be considered is to 
minimize the total system mass while limiting the response to a 
disturbance force at the tip of the Mini-MAST truss. This will be 
accomplished by the selection of the stiffness of the truss elements, the 
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use of a proof mass actuator and a torque wheel actuator, both located at 
the tip. The total system mass is increased when (1) the stiffness of truss 
elements is increased, (2) the capacity of the proof mass actuator is 
increased, and (3) the capacity of the torque wheel actuator is increased. 
Prior to involving the software PDEMOD to calculate the dynamic 
response of the actively controlled Mini-MAST to the disturbance force, 
it is prudent to investigate the best mix of structural stiffness and the use 
of active control in a more simple way. 

The structural stiffness of a uniform Bernoulli beam is given by: 

K = 3EI/L3 

Because the mass of the truss elements is proportional to their stiffness 
the truss mass is: 

MasStruss = MasStruss,o[EI/EI 0 ][L 0 /L]3 

Note that to keep the same resistance to a disturbance force applied at 
the tip of the Mini-MAST the increased stiffness and corresponding mass 
increases by the length to the third power. 

Because the 66-foot Mini-MAST truss has 18 bays it is possible to alter 
the stiffness of the truss elements as each bay to reduce the total mass 
while achieving the same stiffness or resistance to a disturbance force at 
the tip. The maximum saving that can be achieved appears to be about 17 
percent. 

A moving or proof mass actuator (PMA) can produce a force to oppose a 
constant disturbance force but only for a relatively short period of time 
which depends on the size of the proof mass and its stroke. Neglecting 
the stationary mass the mass for the PMA can be shown to be: 

MasspMA = .5*Force*Time 2 /Stroke 

Note that the mass of the proof mass actuator is proportional to the time 
of application squared but does not depend on the truss length or 
stiffness. 

A rotating mass or torque wheel actuator (TWA) can produce a moment 
which by reaction of the truss structure can oppose a constant 
disturbance force but only for a relatively short period. The mass of the 
TWA depends on the moment, maximum wheel speed and the time of 
application. Neglecting the stationary mass the mass of the TWA can be 
expressed as: 

MassxwA = MomentTime/ (Wheel Speed*Radius 2 ) 
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A moment applied to the tip of the Mini-MAST produces a lateral 
deflection which can be used to oppose a disturbance force. The moment 
per deflection is given by: 


Moment /Deflection = 2 EI/L 2 


The force per deflection is: 

Force/Deflection = 3 EI/L 3 

The mass of the totque wheel actuator (TWA) capable of countering a 
particular force for a specific time becomes. 

MassTWA = .667*Force*Time*Length/(Max Wheel Speed*Radius 2 ) 

Note that the mass of the torque wheel actuator is proportional to both 
the time of application of the disturbance force and the length ol the 
Mini-MAST. 

The problem now is to determine the combination of structural stiffness 
and the sizes of the proof mass and torque wheel actuators which 
minimizes the total mass. 

MassTotai = Masstruss + MasspMA + MassTWA 

The fundamental variables which determine the best mix is the time, T, 
for which the disturbance force is applied and the length, L, of the Mini- 
MAST truss. By examining different combinations of truss length, L and 
the time, T, for which the disturbance force is applied the regions for 
which the truss structure, the proof mass actuator, or the torque whee 
actuator is best in countering the disturbance force can be determined. 
The result of such a study shows that the truss structure is bestfor large 
values of the time for which the disturbance force is applied. The prooi 
mass actuator is best for long trusses. Under only limited conditions is 
the torque wheel actuator best for countering disturbance forces. 


CONCLUDING REMARKS 

Distributed parameter models of structures have important advantages for 
problems involving the active control of flexible structures. This is 
especially true for repetitive lattice structures such as the Mini-MAST 
truss because its dynamics can be accurately represented by only a tew 
parameters. 

Until the advent of software for distributed parameter modeling of 
complex flexible structures it was not practical to model complex 
spacecraft configurations using distributed parameter models. Now the 
software PDEMOD enables modeling complex configurations and can also 
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include the control system dynamics in the same equations. The need to 
resort to model order reduction is eliminated because the closed loop 

stability and system performance can be determined without ignoring anv 
of the modes. 6 y 

This capability enables working integrated control/structures problems. 
An example problem is examined in which structure and active controls 
are used to counter the disturbance force at the tip of the NASA Mini- 
MAST truss. The regions are determined in which structure, proof mass 
actuators, and torque wheel actuators result in the minimum mass 
system. Variations of this tradeoff between control and structure 
considerations are being pursued. 

The development of the PDEMOD software has been underway for about 
one year. The formulation and coding has been completed for modelling 
general three-dimensional configurations. Modal frequencies and mode 
shapes have been generated for the Spacecraft Control Laboratory 
Experiment (SCOLE) configuration and the Mini-MAST truss. Graphics 
for drawing wireframes of the deflected configurations has been the most 
recent addition. Transient response, transfer functions, Timoshenko 
beam option and improved root-finding algorithms will be added during 
the next year. Copies of the source code will be made available to anyone 
interested in modelling new configurations or contributing to the 
software development. 
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INTRODUCTION 


• Difficulties of Current Practice 

• Advantages of Distributed Parameter Modeling 

• Difficulties of Modeling Complex Structures 

• Capabilities of PDEMOD Software 

• Integrated Design Objective Functions 

• Necessary Additions to PDEMOD 

• Example Controls-Structures Design Problem 

• Insights Offerred by Particular Tradeoffs 

• Concluding Remarks 
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Finite Element Modeling 

• Excessive Complexity 

• Parameter Estimation is Difficult 

. Model Order Reduction Required for 

Control Analysis 

Distributed Parameter Modeling 

• Fewer Model Parameters 

• Parameter Estimation Straightforward 

. Closed-Loop Stability Analysis does not 
Require Order Reduction 
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Number of Parameters 


Number of Model Parameters 


Solar Array Example: 



15 10 15 20 25 30 35 40 45 50 

Number of Modes 
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Modelin^SoT^art^ 

Structural 

Name Dynamics Optimization 

DISTEL 

by I’oelaert 

BUNVIS 

l)y A rulers on 

PDEMOD 

by Taylor 

+ — Parameter Estimation, 

Control Synthesis, and 
Structural Design 
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New Modeling Software 


• Distributed Parameter Modeling 

• Flex Beam Elements/Rigid Bodies 

(Bend-x, Bcml-y, Twist, Stretch) 

• Parameter Estimation 

• Structural Damping 

• Control System Dynamics 

• Parametric Studies 
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f Beam Model 


y(D 



T(0) 


The Moments and Forces at (0) in Beam Axes are. 


M X = ElyUy(O) 
My =- EI X U^(0) 
M z = EI^ujj,(0) 




F x = EIyUy(0) 
F y = -EI x uf(0) 
F 7 = EA z u£(0) 



i 

M 

2 

1 


Fx 

M 

M v 

F 

Beam 

Fy 

Beam 

y 




_ M z 


_F Z . 


The sign of the forces and moments at the outboard 
end of the beam are opposite to those of the inboard 
end. 
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Generic Configuration 



• Flexible Beams (Bend, Twist, Elongate) 

• Rigid Bodies (Full Inertia Matrices) 
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Partial Differential Equations 
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Bernoul li-Euler Equations will be used for bending 

mu x + El x u ^ =0 

rntiy + EI y u^= o 

String Equations are used for elongation and 
torsion. 

mu z + EA z u z = 0 

Examination of each equation will establish the 
relationship between the b “ parameters and 
frequency, id. 

From ( mu x + El x u x = 0 } we get: 
o 4 

mm z u x = b z EI z u z 


It follows that 




The forces and moments in body axes are: 
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Equations of Motion 


Force Equations: 


A j Q j ^ j^jQ u"j * ffljlD 2 ^^) ^> e am i ®F t 

' J J J body j 


Moment Equations: 


“ ^ + ^2 2t T be ami P M 


body j 


R T p } 

‘'beam i beam i F 
body j body j 


Constraint Equations: 

T 2 U 2 (0> “ R 2 T 2 U 2 (0) " T l u i (G) - R|T|Uj(°) 

T 2 U 2 (0> = T i u/ "i (0) 
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System Equations 


Body 1 
Body 2 
Body 3 
Constraint 


y- Beam 1 

z- Beam 2-, 

/ \ 

Force(w) 

Moment(co) 

0 

Force (co) 
Moment (co) 

Force (coj 
Moment (co) 

0 

Force (co) 
Moment (co) 

U (L) 

u (0) 

u v (L) 

” V (0) 



Characteristic Equation: 

IA(co)l = 0 
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Mode Shape Functions 


u x = A x sinQ3z) + B x cos(pz) + 

C x sinh(Pz) + D x coshQ3z) 


u 





'12 Coefficients 
form State Vector 


Trigonometric and 
Hyperbolic Functions 
^ of o) 



Similarly for u\ Force, Moment, etc. 
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Control System 


r 


V 


A 


Output feedback control systems can be modeled 
directly using the structural model. 

Collocated Force Control 


^bodyj ^u^i + ^u^i 

Ibis controller would add to the appropriate 
term: 

AA M = 4 (K « +|D,K(l3Q u, 

Collocated Moment Control 


^body^ ^1/ 1 

This controller would add to the appropriate 
term: 

AAj j “ 

Higher order compensation for controller can be 
included as complex (unctions of id. 

Uncollocated controllers would have the same form 
but would have different indices to reflect the 
locations of sensors and actuators. 
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Dimensionality 


Because 1 2 modal parameters are invlolved for each beam, A 
like number of equations are involved in the eigen value 
solution The SCOLE configuration is an example of such a simple 
configuration. 

12 


The Mini-MAST, a cantilevered beam with lumped masses at the 
tip and bay 1 0 will involve 24 equations. 

A hexagonal shape with bodies at each point and at its center 
will involve 1 44 equations. 


1 44 


Clearly, numerical difficulties can be expected to limit the 
complexity of a configuration that can be handled. 
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PDEMOD Schedule 


Euler Beam 
Timoshenko 
Modal Freq. 
Mode Shapes 
Wire Frame 
Rigid Body 
Trans. Func. 
Time History 
Systems I.D. 
Damping 
Controls 
Optimization 
Wittrick Alg'm 
Treetops Link 
Applications 



1/90 


1/91 


1/92 


Year 
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Balakrishnan - U.C.LA. 

Huang/ Shen - O.D.U. 

Junkins - Texas A. & M. 

Kakad - U. of N.C. Charlotte 
Kaufman - R.P.I. 

Mitchell - V.P.I. & S.U. 

Pitkin - Dynacs Engineering 
Williams - U. of Nevada @ Las Vegas 
Zi mm erman - U. of Florida 
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Configurations Modeled 




.Spacecraft Control 
Laboratory Experiment 


Mini-MAST— 


Sketches of the Three Spacecraft-Type Structures for which 
Distributed Parameter Models are Constructed. 
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Bending 


Torsion 



Figure 14. Schematics of Distributed Parameter Models for Bending and 
Torsion of the Mini-MAST Truss. 
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Optimization* 


J = 5(q) t m' 1 5(0) 

M = 2 #(e) T M 1 5(0) + 2 1 ^5(0) AO 



* - Parameter Estimation, 
Control Synthesis, 

Structural Design 
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Genetic algorithms (GA's) are one form of directed random search. In GA s a 
finite number of candidate solutions are randomly genereated to create an initial 
population. This initial population is then allowed to evlove overt generations to 
produce new, and hopefully better designs. The basic motivation behind the 
development of GA’s is that they are robust problem solvers for a wide class of 
problems. The basic conjecture behind GA’s is that evolution is the best compromise 
between determinism and chance. 


GENETIC OPTIMIZATION 

o Darwins "Survival of the Fittest" 

o Related to Simulated Annealing 

o Requires Function Evaluations (no Gradient 
Information required) 

o Seeks to Maximize a "Fitness" Index - related to 
objective function 

o Works with Multiple Designs Simultaneously 
o Identifies "Nearly Optimal" Alternatives 
o Suitable for Parallel Processing 
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The genetic algorithm was used to determine the optimal continuum beam 
characteristics of the two sections of the NASA LaRC Mini-Mast for bending in one 
plane of motion. The objective of the design was to minimize the total structural mass 
subject to constraints on the tip displacement. As would be expected, the optimizer 
increases the stiffness of the lower section of the Mini-Mast while decreasing the 
stiffness of the upper section. 


DESIGN PROBLEM 


o FEM Representation of Continuum Model of Mini-Mast 
(PDEMOD soon to follow) 

o Design Variables - EI(z) of Two Sections 
( d*EI(z) nominal ) 

o Problem Statement - Minimize Total Mass Subject to 
Dynamic Displacement Constraints 

o Results 

o Nominal - d(l) = d(2) = 1 

Total Mass = 7.1208 slugs 
Constraint Violation = 5.3% 

o Optimized Design - d(l) = 1.2059 d(2) = 0.4647 

Total Mass = 6.2412 slugs 
Constraint Violation = None 
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The GA search was run for 40 generations with a population size of 20 members. 
Thus, the total nuber of function evaluations was 800. The convergence history o 
GA is shown below At a given generation number, the maximum fitness value 
represents^ The mosT fit member in .be population whereas the average ..ness is the 
Xnfimels of the entire population. I. can be seen that the average , fitness mcreases 
with each new generation, which is a property of the particular GA used. It is also 
seen that the maximum fitness (i.e. optimal solution) is obtained in an early 
generation. 


CONVERGENCE HISTORY 



o Each Generation Represents 20 Function Evaluations 
(corresponding to 20 Designs) 


o Average "Fitness" 
Quickly 


of the Multiple Designs Increase 
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Example Design Problem 

Article: Mini-MAST Truss 

Disturbance: F = 

Specification : u tip < u max 


Controls: 


Objective: 


Proof Mass Actuator 
Torque Wheel Actuator 

J — minfm^russ + m proof 

+ m torque» 
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in Feet 


Control/Structures Optimization 
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Proof Mass Actuator yields lowest Mass 
Torque Wheel Actuator yields lowest Mass 
Structure yields lowest Mass 
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